Studying stellar binary systems with the Laser Interferometer Space Antenna using Delayed Rejection Markov chain Monte Carlo methods

Bayesian analysis of LISA data sets based on Markov chain Monte Carlo methods has been shown to be a challenging problem, in part due to the complicated structure of the likelihood function consisting of several isolated local maxima that dramatically reduces the efficiency of the sampling techniques. Here we introduce a new fully Markovian algorithm, a Delayed Rejection Metropolis-Hastings Markov chain Monte Carlo method, to efficiently explore these kind of structures and we demonstrate its performance on selected LISA data sets containing a known number of stellar-mass binary signals embedded in Gaussian stationary noise.Comment: 12 pages, 4 figures, accepted in CQG (GWDAW-13 proceedings

A Bayesian parameter estimation approach to pulsar time-of-arrival analysis

The increasing sensitivities of pulsar timing arrays to ultra-low frequency (nHz) gravitational waves promises to achieve direct gravitational wave detection within the next 5-10 years. While there are many parallel efforts being made in the improvement of telescope sensitivity, the detection of stable millisecond pulsars and the improvement of the timing software, there are reasons to believe that the methods used to accurately determine the time-of-arrival (TOA) of pulses from radio pulsars can be improved upon. More specifically, the determination of the uncertainties on these TOAs, which strongly affect the ability to detect GWs through pulsar timing, may be unreliable. We propose two Bayesian methods for the generation of pulsar TOAs starting from pulsar "search-mode" data and pre-folded data. These methods are applied to simulated toy-model examples and in this initial work we focus on the issue of uncertainties in the folding period. The final results of our analysis are expressed in the form of posterior probability distributions on the signal parameters (including the TOA) from a single observation.Comment: 16 pages, 4 figure

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.Comment: 12 pages, 2 figures, proceedings of the 8th Edoardo Amaldi Conference on Gravitational Waves, New York, June 21-26, 200

The Mock LISA Data Challenges: from Challenge 1B to Challenge 3

The Mock LISA Data Challenges are a programme to demonstrate and encourage the development of LISA data-analysis capabilities, tools and techniques. At the time of this workshop, three rounds of challenges had been completed, and the next was about to start. In this article we provide a critical analysis of entries to the latest completed round, Challenge 1B. The entries confirm the consolidation of a range of data-analysis techniques for Galactic and massive--black-hole binaries, and they include the first convincing examples of detection and parameter estimation of extreme--mass-ratio inspiral sources. In this article we also introduce the next round, Challenge 3. Its data sets feature more realistic waveform models (e.g., Galactic binaries may now chirp, and massive--black-hole binaries may precess due to spin interactions), as well as new source classes (bursts from cosmic strings, isotropic stochastic backgrounds) and more complicated nonsymmetric instrument noise.Comment: 20 pages, 3 EPS figures. Proceedings of the 12th Gravitational Wave Data Analysis Workshop, Cambridge MA, 13--16 December 2007. Typos correcte

First LIGO search for gravitational wave bursts from cosmic (super)strings

We report on a matched-filter search for gravitational wave bursts from cosmic string cusps using LIGO data from the fourth science run (S4) which took place in February and March 2005. No gravitational waves were detected in 14.9 days of data from times when all three LIGO detectors were operating. We interpret the result in terms of a frequentist upper limit on the rate of gravitational wave bursts and use the limits on the rate to constrain the parameter space (string tension, reconnection probability, and loop sizes) of cosmic string models.Comment: 11 pages, 3 figures. Replaced with version submitted to PR

A Joint Search for Gravitational Wave Bursts with AURIGA and LIGO

The first simultaneous operation of the AURIGA detector and the LIGO observatory was an opportunity to explore real data, joint analysis methods between two very different types of gravitational wave detectors: resonant bars and interferometers. This paper describes a coincident gravitational wave burst search, where data from the LIGO interferometers are cross-correlated at the time of AURIGA candidate events to identify coherent transients. The analysis pipeline is tuned with two thresholds, on the signal-to-noise ratio of AURIGA candidate events and on the significance of the cross-correlation test in LIGO. The false alarm rate is estimated by introducing time shifts between data sets and the network detection efficiency is measured with simulated signals with power in the narrower AURIGA band. In the absence of a detection, we discuss how to set an upper limit on the rate of gravitational waves and to interpret it according to different source models. Due to the short amount of analyzed data and to the high rate of non-Gaussian transients in the detectors noise at the time, the relevance of this study is methodological: this was the first joint search for gravitational wave bursts among detectors with such different spectral sensitivity and the first opportunity for the resonant and interferometric communities to unify languages and techniques in the pursuit of their common goal.Comment: 18 pages, IOP, 12 EPS figure

All-sky search for periodic gravitational waves in LIGO S4 data

We report on an all-sky search with the LIGO detectors for periodic gravitational waves in the frequency range 50-1000 Hz and with the frequency's time derivative in the range -1.0E-8 Hz/s to zero. Data from the fourth LIGO science run (S4) have been used in this search. Three different semi-coherent methods of transforming and summing strain power from Short Fourier Transforms (SFTs) of the calibrated data have been used. The first, known as "StackSlide", averages normalized power from each SFT. A "weighted Hough" scheme is also developed and used, and which also allows for a multi-interferometer search. The third method, known as "PowerFlux", is a variant of the StackSlide method in which the power is weighted before summing. In both the weighted Hough and PowerFlux methods, the weights are chosen according to the noise and detector antenna-pattern to maximize the signal-to-noise ratio. The respective advantages and disadvantages of these methods are discussed. Observing no evidence of periodic gravitational radiation, we report upper limits; we interpret these as limits on this radiation from isolated rotating neutron stars. The best population-based upper limit with 95% confidence on the gravitational-wave strain amplitude, found for simulated sources distributed isotropically across the sky and with isotropically distributed spin-axes, is 4.28E-24 (near 140 Hz). Strict upper limits are also obtained for small patches on the sky for best-case and worst-case inclinations of the spin axes.Comment: 39 pages, 41 figures An error was found in the computation of the C parameter defined in equation 44 which led to its overestimate by 2^(1/4). The correct values for the multi-interferometer, H1 and L1 analyses are 9.2, 9.7, and 9.3, respectively. Figure 32 has been updated accordingly. None of the upper limits presented in the paper were affecte

Search for gravitational waves from binary inspirals in S3 and S4 LIGO data

We report on a search for gravitational waves from the coalescence of compact binaries during the third and fourth LIGO science runs. The search focused on gravitational waves generated during the inspiral phase of the binary evolution. In our analysis, we considered three categories of compact binary systems, ordered by mass: (i) primordial black hole binaries with masses in the range 0.35 M(sun) < m1, m2 < 1.0 M(sun), (ii) binary neutron stars with masses in the range 1.0 M(sun) < m1, m2 < 3.0 M(sun), and (iii) binary black holes with masses in the range 3.0 M(sun)< m1, m2 < m_(max) with the additional constraint m1+ m2 < m_(max), where m_(max) was set to 40.0 M(sun) and 80.0 M(sun) in the third and fourth science runs, respectively. Although the detectors could probe to distances as far as tens of Mpc, no gravitational-wave signals were identified in the 1364 hours of data we analyzed. Assuming a binary population with a Gaussian distribution around 0.75-0.75 M(sun), 1.4-1.4 M(sun), and 5.0-5.0 M(sun), we derived 90%-confidence upper limit rates of 4.9 yr^(-1) L10^(-1) for primordial black hole binaries, 1.2 yr^(-1) L10^(-1) for binary neutron stars, and 0.5 yr^(-1) L10^(-1) for stellar mass binary black holes, where L10 is 10^(10) times the blue light luminosity of the Sun.Comment: 12 pages, 11 figure

All-sky LIGO Search for Periodic Gravitational Waves in the Early S5 Data

We report on an all-sky search with the LIGO detectors for periodic gravitational waves in the frequency range 50--1100 Hz and with the frequency's time derivative in the range -5.0E-9 Hz/s to zero. Data from the first eight months of the fifth LIGO science run (S5) have been used in this search, which is based on a semi-coherent method (PowerFlux) of summing strain power. Observing no evidence of periodic gravitational radiation, we report 95% confidence-level upper limits on radiation emitted by any unknown isolated rotating neutron stars within the search range. Strain limits below 1.E-24 are obtained over a 200-Hz band, and the sensitivity improvement over previous searches increases the spatial volume sampled by an average factor of about 100 over the entire search band. For a neutron star with nominal equatorial ellipticity of 1.0E-6, the search is sensitive to distances as great as 500 pc--a range that could encompass many undiscovered neutron stars, albeit only a tiny fraction of which would likely be rotating fast enough to be accessible to LIGO. This ellipticity is at the upper range thought to be sustainable by conventional neutron stars and well below the maximum sustainable by a strange quark star.Comment: 6 pages, 1 figur

Search for Gravitational Waves Associated with 39 Gamma-Ray Bursts Using Data from the Second, Third, and Fourth LIGO Runs

We present the results of a search for short-duration gravitational-wave bursts associated with 39 gamma-ray bursts (GRBs) detected by gamma-ray satellite experiments during LIGO's S2, S3, and S4 science runs. The search involves calculating the crosscorrelation between two interferometer data streams surrounding the GRB trigger time. We search for associated gravitational radiation from single GRBs, and also apply statistical tests to search for a gravitational-wave signature associated with the whole sample. For the sample examined, we find no evidence for the association of gravitational radiation with GRBs, either on a single-GRB basis or on a statistical basis. Simulating gravitational-wave bursts with sine-gaussian waveforms, we set upper limits on the root-sum-square of the gravitational-wave strain amplitude of such waveforms at the times of the GRB triggers. We also demonstrate how a sample of several GRBs can be used collectively to set constraints on population models. The small number of GRBs and the significant change in sensitivity of the detectors over the three runs, however, limits the usefulness of a population study for the S2, S3, and S4 runs. Finally, we discuss prospects for the search sensitivity for the ongoing S5 run, and beyond for the next generation of detectors.Comment: 24 pages, 10 figures, 14 tables; minor changes to text and Fig. 2; accepted by Phys. Rev.